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The article reviews recent progress in the theoretical understanding of near-field surface electro- 
magnetic phenomena in pristine and atomically doped carbon nanotubes. The phenomena involv- 
ing strong coupling effects are outlined. They are the optical absorption by single-walled carbon 
nanotubes doped with single atoms or ions in the frequency range close to the atomic transition 
frequency, the entanglement of the pair of atomic qubits strongly coupled to a common high-finesse 
surface photonic mode of the nanotube, and the optical response of the strongly coupled surface 
exciton-plasmon excitations in pristine semiconducting carbon nanotubes. The phenomena reviewed 
have a great potential to be exploited for the future development of the nanotube based tunable 
optoelectronic device applications in areas such as nanophotonics, nanoplasmonics, cavity quantum 
electrodynamics, and quantum information science. 
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1. Introduction 

Single-walled carbon nanotubes (CNs) are quasi-one- 
dimensional (ID) cylindrical wires consisting of graphcnc 
sheets rolled- up into cylinders with diameters ~ 1 — 10 nm 
and lengths - 1 - 10 4 /mi Q, H H [|. Nanotubes 
are shown to be useful as miniaturized electronic, elec- 
tromechanical, and chemical devices [f|, scanning probe 
devices 0, and nanomaterials for macroscopic compos- 
ites Q- Their intrinsic physical properties may be 
substantially modified in a controllable way by dop- 
ing with extrinsic impurity atoms, molecules and com- 
pounds @, 0, M, M\mM\u\- 

Recent successful experiments on the encapsulation of 
single atoms into single- walled CNs [IH, Il3l.[l3 and their 
intercalation into single- wall CN bundles [lfj. IH]. along 
with the numerous previous studies of monoatomic gas 
absorption by the nanotube bundles (see, e.g., Ref. [151 ]) 
and the progress in growth techniques of centimeter-long 
small-diameter single- walled carbon nanotubes [H, 
stimulated an in-depth theoretical analysis of dynamic 
quantum coherent processes in CNs doped with single 
atoms or ions. As a result, atomically doped nanotubes 
have been recently demonstrated to be very efficient as 
applied to nanophotonics, cavity quantum electrodynam- 
ics (QED) and quantum communication 1|| 17, Tg> 19j - 
For pristine (undoped) single-walled CNs, the numer- 
ical calculations predicting large exciton binding ener- 
gie s (~ 0.3 — 0.6 eV) in semiconducting CNs [2CJ, [2l|, 
|22| and even in some small-diameter (~ 0.5 nm) metal- 
lic CNs 0, followed by the results of various exciton 



photoluminescence measurements 0, H^, 2|| 2?], [H, [2i| , 
have become available. These works, together with other 
reports investigating the role of effects such as intrinsic 
defects [26J, |30[ , exciton-phonon interactions [27L m hi 
32, 33j], exciton-surface-plasmon coupling 34, 39 ] , eff ects 
of external magnetic [2g, H3] and electric fields [38| , re- 
veal the variety and complexity of the intrinsic optical 



properties of semiconductor carbon nanotubes [391 ]. At 
the same time, they pave the way for the novel applica- 
tions of carbon nanotubes in modern optoelectronics in 
areas such as various aspects of photophysics, nanoplas- 
monics, and quantum information science, which are es- 
sentially beyond relatively well explored nanoelectronics 
and nanoelectromechanics with carbon nanotubes. For 
instance, recent experimental observations of quantum 
correlations in the photoluminescence spectra of individ- 
ual CNs suggest their potential applicability in quantum 
cryptography [29j |. 

In spite of impressive experimental demonstrations of 
basic quantum information effects in a number of dif- 
ferent nanoscale solid state systems, such as quantum 
dots in semiconductor microcavities, nuclear spin sys- 
tems, Josephson junctions, etc., their concrete implemen- 
tation is still at the proof-of-principle stage (see, e.g., 
Ref. [i^] for a review). The development of materials 
that may host quantum coherent states with long co- 
herence lifetimes is a critical research problem for the 
nearest future. There is a need for the fabrication of 
quantum bits (qubits) with coherence lifetimes at least 
three-four orders of magnitude longer than it takes to 
perform a bit flip [i^, EH, To achieve this goal, a 
critical prerequisite of the strong coupling of an atomic, 
excitonic or intersubband electronic transition to a high- 
finesse vacuum-type electromagnetic mode of its environ- 
ment (e.g., a semiconductor optical microcavity mode) 
needs to be fulfilled. Under strong coupling, new ele- 
mentary excitations, the eigen states of the full photon- 
matter Hamiltonian, form. They are so-called 'cavity po- 
laritons' representing a characteristic anticrossing behav- 
ior with a mode separation often referred to as 'vacuum- 
field Rabi splitting'. A pair of the cavity polaritons 
(term used in a broad sense here) coupled to the same 
high-finesse vacuum-type electromagnetic mode of the 
surrounding material would involve polariton entangling 
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operations 17|, |4l|, |42|, |43j, |44|, followed by the nearest 
neighbor interaction over short distances and quantum 
information transfer over longer distances. 

For atomically doped carbon nanotubes, there may be 
typically two qualitatively different regimes of the in- 
teraction of an atomic state with the medium-assisted 
surface electromagnetic modes of a CN j4|| 46, 47, 48, 
l49l . IHoj . They are the weak coupling regime and the 
strong coupling regime. The former is characterized by 
the monotonic exponential decay dynamics of the up- 
per atomic state with the decay rate altered compared 
to the free-space value due to the presence of the ad- 
ditional decay modes associated with the presence of 
the nanotube. The latter is, in contrast, strictly non- 
exponential and is characterized by the reversible (Rabi) 
oscillations where the energy of the initially excited atom 
is periodically exchanged between the atom and the sur- 
face electromagnetic mode of the nanotube. The strong 
coupling regime was shown recently to occur in small- 
diameter (~ 1 nm) atomically doped CNs with the atomic 
transition tuned to the resonance (representing the field 
mode) of the local density of the nanotube's surface pho- 



tonic states (DOS) [46j, |47J, |4& |49|, |5C(. This result 



can be qualitatively understood in terms of the cavity 
QED. The coupling constant of an atom (modelled by 
a two-level system with the transition dipole moment 
d,A and frequency loa) to a local vacuum- type field is 
given by fig — (2Trd' A hujA/V) 1 ^ 2 with V being the effec- 
tive volume of the field mode the atom interacts with 
(see, e.g., Ref. [Hl[). For the atom (ion) encapsulated 
into the CN of radius R cn , the local field mode volume 
is estimated to be V ~ itR^^Xa/Z) that is ~ 10 2 nm 3 
for CNs with diameters ~ 1 nm in the optical ran ge o f 
Xa ~ 600 nm. Approximating d,A ~ er ~ e(e 2 /hu)A) [52| . 
one obtains %~0.3 eV. On the other hand, the "cavity" 
linewidth is given for u>a in resonance with the cavity 
mode by hjc = 6irhc 3 /ui A t;(ujA)V , where £ is the atomic 
spontaneous emission enhancement/dehancement (Pur- 
cell) factor [Hl|. Takin g in to account large Purcell fac- 
tors ~10 7 close to CNs [45[, one arrives at ?i7 c ^0.03 eV 
for 1 nm-diameter CNs in the optical spectral range. 
Thus, for the atoms (ions) encapsulated into the small- 
diameter CNs the strong atom-field coupling condition 
9 lie 3> 1 is fulfilled, giving rise to the rearrangement 
("dressing") of the atomic levels followed by the forma- 
tion of q uasi-one-dimensional (ID) atomic cavity polari- 



tons [S E3, IS [H3l. The latter ones are similar to 
quasi-OD excitonic polaritons in quantum dots in semi- 

[Ei, " 



conductor microcavities 



which arc theoreti- 



cally proposed to be a way to produce the excitonic states 
entanglement for use in solid state quantum information 
science applications |43j |. Similar way to entangle the 
pair of spatially separated quasi- ID atomic polaritons in 
metallic small-diameter single-walled CNs has been the- 
oretically demonstrated to be feasible as well 17]. It is 
important to pursue a variety of different strategies and 



approaches towards physically implementing novel non- 
trivial applications in modern nanotechnology. 

For pristine small-diameter (< 1 nm) semiconducting 
single-walled CNs, the exciton-surface-plasmon interac- 
tions have been shown recently to result in the strong 
exciton-plasmon coupling regime [H, [35| . This is due to 
the fact that the low-energy weakly-dispersive interband 
plasmon modes (resulted from the transverse quantiza- 
tion of the electronic motion on the nanotube surface — 
observed in Ref. [5^]) occur in the same energy range of 
~ 1 eV where the lowest exciton excitation energies are 
located in small-diameter CNs [EH. Previous stud- 
ies of the exciton-plasmon interactions in nanomaterials 
have been focused on artificially fabricated hybrid plas- 
monic nanostructures, such as dye molecules in organic 
polymers deposited on metallic films [59j |. semiconduc- 
tor quantum dots coupled to metallic nanoparticles [ooj ] . 
or nanowires 6l| , where one material carries the exciton 
and another one carries the plasmon. The effect of the 
strong exciton-plasmon coupling in CNs is particularly 
interesting since here the fundamental electromagnetic 
phenomenon of strong coupling occurs in an individual 
quasi-lD nanostructure, a carbon nanotube. 

This article briefly reviews recent progress in the the- 
oretical understanding of near-field surface electromag- 
netic phenomena in pristine and atomically doped CNs. 
It focuses on the phenomena that involve the strong 
coupling effects. Specifically, light absorption by atomi- 
cally doped CNs and the entanglement of spatially sep- 
arated strongly coupled quasi-lD atomic polaritons are 
discussed. Also, physical reasons for the exciton-surface- 
plasmon coupling in pristine small-diameter semicon- 
ducting CNs is given, and the optical response of the 
strongly coupled surface exciton-plasmon excitations is 
analyzed. These strong coupling phenomena have a great 
potential to be exploited for the future development of 
the nanotube based tunable optoelectronic device appli- 
cations in areas such as nanophotonics, nanoplasmonics, 
cavity QED, and quantum information science. 

The macroscopic electromagnetic Green function for- 
malism is used throu gho ut the paper, developed in 
Refs. jlil l4tjL l47l lij. l49l i50ll to study the quantum nature 
of the near-field electromagnetic interactions in quasi- ID 
nanostructures featuring strong absorption. The formal- 
ism follows the original line of the macroscopic QED ap- 
proach introduced by Welsch and coworkers to rigorously 
describe medium- assisted electromagnetic vacuum effects 
in dispersing and absorbing media 62, 63, 6^] (and refs. 
therein) . 

The paper is organized as follows. Section 2 ana- 
lyzes the optical absorption by single-walled CNs doped 
with single atoms (or ions) in the frequency range close 
to the atomic transition frequency. A significant line 
(Rabi) splitting effect is demonstrated for small-diameter 
(~ 1 nm) CNs 16[ , which indicates the strong atom- field 
coupling regime known to facilitate the entanglement of 
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spatially separated atomic qubits 4J, |65| . Section 3 cal- 
culates the entanglement of formation for the pair of 
atomic qubits (quasi-lD atomic polaritons) strongly cou- 
pled to the same medium-assisted high-finesse surface 
photon mode of the nanotube. Small-diameter metal- 
lic nanotubes are shown to result in sizable amounts of 
the two-qubit atomic entanglement for sufficiently long 
times [17| . Section 4 describes the physical nature of the 
interactions between excitonic states and interband sur- 
face plasmons in pristine small-diameter semiconducting 
single-walled CNs, and calculates the optical response 
(exciton absorption) under strong exciton-plasmon cou- 
pling. Section 5 discusses the envisaged nanotube based 
applications of the phenomena considered, summarizes 
and concludes the article. 

2. Optical Absorbtion by Atomically Doped CNs 

The (achiral) CN is modelled by an infinitely long, in- 
finitely thin, anisotropically conducting cylinder. Its (ax- 
ial) surface conductivity is taken to be that given by the 
7r-electron band structure in the tight-binding approx- 
imation with the azimuthal electron momentum quan- 
tization and axial electron momentum relaxation taken 
into account. A two- level atom is positioned at the point 
ta near an infinitely long achiral single- wall CN (Fig.[T]). 
The orthonormal cylindric basis {e r , e v , e z } is chosen in w [th 
such a way that e z is directed along the nanotube axis 
and, without loss of generality, ya = ^46 r — {^,0,0}. 
The atom interacts with the quantum electromagnetic 
field via its transition dipole moment directed along the 
CN axis, dA =d z e z . The contribution of the transverse 
dipole moment orientations is suppressed because of the 
strong depolarization of the transverse field in an isolated 



CN [66J, [67|, |68|, 16J, HO, UK]. Under these assumptions, 



the photon Green function formalism developed earlier 
in Refs. 0, for quantizing an electromagnetic field 
close to quasi- ID absorbing and dispersing media gener- 
ates the total second quantized Hamiltonian (Gaussian 
units) 

H = J du> huJdR /t(R, w)/(R, u) + ^ a z (1) 

+ J^dujJdR [gW^^.w)^ 

-g(-)(r^,R,o;)a]/(R,a;)+h.c., 

where the three terms represent the medium-assisted 
(modified by the presence of the CN) electromagnetic 
field, the two-level atom and their interaction, respec- 
tively. The operators P(R, u>) and /(R, u>) are the scalar 
bosonic field operators defined on the CN surface as- 
signed by the radius- vector R = {i? crl , 0, Z} with R cn be- 
ing the radius of the CN. These operators create and an- 
nihilate the single-quantum bosonic-type electromagnetic 



medium excitation of the frequency u> at the point R of 
the CN surface. The Pauli operators, a z =\u){u\ — 
a= \l){u\ and = \u)(l\, describe the atomic subsystem 
and electric dipole transitions between the two atomic 
states, upper \u) and lower \l), separated by the transi- 
tion frequency lja- This (bare) frequency is modified by 
the diamagnetic (~ A 2 ) atom- field interaction yielding 
the new renormalized frequency 



UJA = UJA 



1 



J^dujJdR^irA^ujtf 



(hu A ) 2 . 

in the second term of the Hamiltonian. The dipole atom- 
field interaction matrix elements g^^r^, R, oj) are given 
by g ( &=g L ± {ui/uj A )gK where 

-(id. 



x ^irh,ojRe(T zz (oj) ""G zz (rA, R, w) 

with ± ^G ZZ being the zz-component of the transverse 
(longitudinal) Green tensor (with respect to the first vari- 
able) of the electromagnetic subsystem, and a zz (u>) rep- 
resenting the CN surface axial conductivity. The matrix 
elements g^dD have the property of 

t. 

2n V oj ) 



dR|g- 



-(||>(^,R, W )| 2 = M^) r H£- 



.(II) 



^ {ll) (r A ,uj) 



Im^G^irA^A^) 



ImG°» 

being the transverse (longitudinal) distance-dependent 
(local) photonic DOS functions, and 

8TTUJ 2 d 2 



T (oj) 



3hc 2 



ImG°» 



representing the atomic spontaneous decay rate in vac- 
uum, where \to.G° zz {uj)=uj /6ttc is the vacuum imaginary 
Green tensor zz-component. 

The Hamiltonian |T]) involves only two standard ap- 
proximations. They are the electric dipole approxima- 
tion and the two-level approximation. The rotating wave 
approximation commonly used is not applied, and the 
diamagnetic term of the atom-field interaction is not ne- 
glected (as opposed to, e.g., Refs. 4^, 46, 47, 65|). Quan- 
tum electrodynamics of the two-level atom (ion) close to 
the CN is thus described in terms of only two intrin- 
sic characteristics of the electromagnetic subsystem - 
the transverse and longitudinal local photonic DOS func- 
tions. 

When the atom is initially in the upper state and the 
field subsystem is in vacuum, the time-dependent wave 
function of the whole system can be written as 

\m) =Cu(t)e-^^ t \u)\{0}) 



(2) 



rfR / dujCi(R,u),t)e 



— i(uj— u)£_/'2)t I 



|l(R,w)>, 
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FIG. 1: (Color online) The geometry of the problem. 

where |{0}) is the vacuum state of the field subsystem, 
|{1(R, cl>)}) is its excited state with the field being in the 
single-quantum Fock state, C u and C\ are the population 
probability amplitudes of the upper state and the lower 
state of the whole system, respectively. Applying the 
Hamiltonian |T]) to the wave function ([2]), one has 

C u (t) = --jdu;jdTlg( + Hr A ,Il,u J ) (3) 



Q (R, w, t) = - ± [g<+> (r A , R, «)r C u (t) 



i(u)—u>A)t 



(4) 



In terms of the probability amplitudes above, the 
emission/ absorbtion spectral line shape I{u>) is given by 
JdR\Ci(H,oj,t -> oo)| 2 [z3. (Obviously, the absorb- 



tion line shape coincides with the emission line shape if 
the monochromatic incident light beam is used in the 
absorbtion experiment.) This, after the substitution of 
Ci obtained by the integration [under initial conditions 
Cj(R,a;,0) = 0, Cu(0) = 1] of Eq. |@J), yields in dimen- 
sionless variables 



I(r A ,x) = I (r A ,x A ) 



where 



2 7o 7 



r = 



270 



drCu(T) e«' 



TlLJ 

270" 



(5) 



2 7o i 



with 7o = 2.7 eV being the carbon nearest neighbor hop- 
ping integral (ft/270 = 1.22 x 10 -16 s) appearing in the 
CN surface axial conductivity a zz , and 



Io(r A ,x) 



f Qr) 
2tt 



(^) 2 ^(r A , X )+tHr A ,x) 



V X ) 



The upper state population probability amplitude in 
Eq. ([5]) is given by the Volterra integral equation [ob- 
tained by substituting the result of the formal integration 



of Eq. ([4]) into Eq. ([3])] with the kernel determined by the 
local photonic DOS functions ^ n (r A , x) 0, 0, [H^. 
Thus, the numerical solution is only possible for the line 
shape /, strictly speaking. This, however, offers very 
little physical insight into the problem of the optical ab- 
sorbtion by atomically doped CNs under different atom- 
field coupling regimes. On the other hand, for those 
atomic transition frequencies x A which are located in the 
vicinity of the resonance frequencies x r of the DOS func- 
tions £-"-(11) a simple analytical approach may be applied. 
Specifically, £,^^{y a , x~x r ) can be approximated by the 
Lorentzians of the same half- width- at-half-maxima 6x r , 
thus making it possible to solve the integral equation for 
C u analytically to obtain 



Sx r 



V / Sx 2 -X 2 / 

Sx r 
^/5x 2 r - X 2 



(6) 



with 



X 



\J 25x r r (x A )^(r A ,x A ) . 



(7) 



This solution is valid for x A f=a x r whatever the atom- 
field coupling strength is, yielding the exponential decay 
of the upper state population, |Cw(t)| 2 ~ exp[— T(x a )t] 
with the rate T f=a Fo^ -1 , in the weak coupling regime 
where (X/5x r ) 2 <§C 1, and the decay via damped Rabi 
oscillations, \C u (t)\ 2 ~ exp(— Sx r r) cos 2 (Xt/2), in the 
strong coupling regime where (X/Sx r ) 2 3> 1. 

Substituting Eq. ([6]) into Eq. (J5]) and integrating 
over t, one arrives at the expression 



I (r A , x A ) [(x - x A ) 2 + 5x1] 
[(x - x A ) 2 - X 2 /A} 2 + Sx 2 {x - x A f 



(8) 



representing the emission/absorbtion spectral line shape 
for frequencies x ~ x r ~ x A regardless of the atom-field 
coupling strength. The line shape ((8]) is clearly seen to be 
of a symmetric two-peak structure in the strong coupling 
regime where (X/Sx r ) 2 ^> 1. The exact peak positions 
are 



£1,2 = x A ± 



(A72)/yr 



+ 8{Sx r /X) 2 - A{5x r /Xf 



separated from each other by x\ — X2 ~ X with the 
Rabi frequency X given by Eq. 0. In the weak cou- 
pling regime (X/Sx r ) 2 <C 1, and xi t 2 become complex, 
indicating that there are no longer peaks at these fre- 
quencies. As this takes place, Eq. JHJ) is approximated 
with the weak coupling condition, the fact that x ~ x A 
and Eq. ([7]), to give the well-known Lorentzian 



I(x) 



(x~x A ) 2 +T 2 (x A )/4 
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FIG. 2: (Color online) Dimensionless (in units of e 2 /2irTi) 
axial surface conductivity of the (11,0) CN per unit length. 
Dimensionless energy is defined as [Energy] /2 r yo- 



with the half-width-at-half-maximum T(xa)/2, peaked 
at x —— x a • 

To compute the absorption line shapes of particular 
atomically doped CNs, one needs to know their trans- 
verse and longitudinal local photonic DOS functions 
i^\Y Al x A ) in Eqs. ® and Q. These arc determined 
by the CN surface axial conductivity a zz , which is to be 
calculated beforehand. Figure [5] shows an example of the 
surface axial conductivity for the (11,0) nanotube 13511 . 
calculated in the random-phase approximation [zHTzil 
with the CN band structure taken into account within 
the nearest-neighbor non-orthogonal tight-binding ap- 
proach [75l ]. 

Figure [3] (a) shows the functions £ (x) for the atom 
in the center of the (5,5), (10,10) and (23,0) CNs (calcu- 



lated as described in Refs. 0, Sfl |47|, |4& M, |50j ) . They 
are seen to decrease with increasing the CN radius, repre- 
senting the decrease of the atom-field coupling strength 
as the atom moves away from the CN wall. To calcu- 
late the normalized absorption curves A(x) = I (x) / I pea k 
in this case, the atomic transition frequency was fixed 
at xa — 0.45 [vertical dashed line in Fig. [3] (a)] since 
this is the approximate peak position of £ for all the 
three CNs. Figure [3] (b) shows the result. For the 
(5,5) CN, the absorption line is split into two compo- 
nents, indicating the strong atom-field coupling with the 
Rabi splitting X « 0.2 corresponding to the energy of 
0.2 x 270 = 1.08 eV. The splitting is disappearing with in- 
creasing the CN radius, and the absorption line shape be- 
comes Lorentzian with the narrower widths for the larger 
radii CNs. 

3. Atomic States Entanglement in CNs 



Two identical two- level atoms (pair of atomic qubits), 
A and B, located at va ={tai 0, 0} and vb = {tb, <Pb, zb} 
close to a CN, are coupled to the common medium- 
assisted vacuum- type electromagnetic field via their tran- 
sition dipole moments. Both atomic dipole moments are, 
as before, assumed to be directed along the CN axis, 
= ds = d z e z . Also assumed is that the atoms are lo- 
cated far enough from each other, to simplify the problem 
by ignoring the interatomic Coulomb interaction. The to- 
tal secondly quantized Hamiltonian of the system of the 
two identical atoms and the nanotube is then written as 
follows [I3| [compare with Eq. |T|)] 



H = 



Tiujj 



dw huo / dR p (R, w)/(R, lu) + J2 -7; ~ *i» (9) 



i=A,B 



+ V / < L/dR[gW(r <> R,a;)o-t 

i=A,B J ° J 

-g ( - ) (r J ,R,w)a l ]/(R,^)+h.c.. 

The Pauli operators, <T; = a\ — \ui)(L\ and &i Z = 

\v>i)(ui \ — [L)(L[ with i = A or B, describe the atomic 




r A =0, x A =0.45 

-■- (5,5) 

(10,10) 

(23,0) 



FIG. 3: (Color online) Transverse local photonic DOS (a) 
and normalized absorption line shapes (b) for the atom in the 
center of the CNs of increasing radii. See text for notations. 
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subsystem and electric dipole transitions between its two 
states, upper \ui) with either of the two atoms in its 
upper state and lower \L) with both atoms in their lower 
states, separated by the transition frequency w A . The 
renormalized transition frequencies u>a,b are the same 
for both of the atoms. 

For single-quantum excitations, the time-dependent 
wave function of the whole system can be written as 

w)) =E^) e " (M)t K->K Q }> ( 10 ) 

i=A,B 

+ dujJdR. C L (R, u, t) e-^-"^ \L) |1(R, u)) , 

where Zu =J2i=A b &i/2 = &A, CVt; and Cl are the proba- 
bility amplitudes for the upper states and the lower state 
of the system, respectively. For the following it is conve- 
nient to introduce the new variables 

C ± {t) = ±=[Cu A {t)±C UB {t)] 

that are the expansion coefficients of the wave func- 
tion (jlOp in terms of the maximally entangled two-qubit 
atomic states \±} = (\u A ) ±\u B ))/\/2. In view of Eqs. ([9|) 
and (fl0|) , the time-dependent Schrodinger equation yields 
then (dimensionless variables — Section 2) 



C±(t) = / cLt'K ± {t - t') C±(t') + /±(r) , (11) 
Jo 



are the generalizations of the DOS functions f-'-dl) (r^, x) 
of Section 2 (see Refs. [13, 0, 0, H3] for more details). 



where 



K±(t-t') = - — / dxT (x) 
27T Jo 



(12) 



x ^ ± {^A^B,x)e 



-i{x — XA){T-T') 



^(r A ,r B ,x) = ^ [^{r Al x)±^{v A ,v Bl x)} (13) 



-^{r A ,x)±^ l {r A ,r B ,x) 



f±{r) = -^=fdT'K ± (T-r')G l 



■u A (r'). (14) 



The functions /±(t) are only unequal to zero when the 
two atoms are initially in their ground states, with the 
initial excitation residing in the nanotube. Equation ([Til) 
assumes that this is realized by selecting the time origin 
to be right after the time interval At, that is necessary 
for the (excited) atom A to decay completely into the 
nanotube photonic modes, has elapsed. The two-particle 
local photonic DOS functions 



e ill) (r A ,r B ,x) 



, ; >/ with ay'" 8 " 1 be- 



The entanglement of the two quantum bits occurs 
when the two-qubit wave function cannot be represented 
as a product of the two one-qubit states in any basis. 
To determine this quantity, a recipe based on the "spin 
flip" transformation (see Ref. (7(| for details) and valid 
for an arbitrary number of the qubits is used. First, the 
reduced density matrix is defined 

Pab(t) = \iP AB (t))(^ ab (t)\ = Tr /ieid |-0(r))(-0(r)| , 

that describes the bipartite atomic subsystem in terms 
of the wave function (fTQ|) of the whole system. Next, the 
" concurrence" 

G(4> AB ) = \ (iPab\iPab}\ 

is introduced, where \iPab) = < J y®y\ > l } 
ing the Pauli matrix that represents the " spin flip" trans- 
formation in the atom A(B) single-qubit space. This, af- 
ter some algebra, becomes 

C^ AB (r)} = \Cl(r)-C 2 _(r)\ (15) 

with C±(r) given by the integral equation (fTTj) . Finally, 
the degree of the entanglement of the two-qubit atomic 
state \ip AB ) ("entanglement of formation" [7y|) is given 
by 

E[^ B (r)] = ft jl + - C[^b(t)] 2 | , (16) 

where h(y) = -y log 2 y - (1 - y) log 2 (l - y). 

The entanglement Ei[tp AB (r)] is maximal when the co- 
efficients C+ (r) and C_ (r) are maximally different from 
each other. For this to occur, the functions £ ( r Aj r B , x) 
in Eqs. (| 1 1 [) — (fl4|) should be different in their values. 
These are determined by £ ± ^ (r A , r B: x) whose fre- 
quency behavior is determined by the CN surface axial 
conductivity a zz (Fig. [5]). The functions ^(H) (r^, r B , x) 
are computed as described in Ref. [13]. Then, the inte- 
gral equation (fTTj) is solved numerically to obtain C±(r), 
followed by the " concurence" (fT5"]) and the entanglement 
of formation (JTHJ) . 

Shown in Fig. U (a) is a typical frequency behavior of 
the one-particle transverse photonic DOS ^(r-A = 0, x) 
in the infrared and visible frequency range x < 0.4 for the 
atom inside 'zigzag' CNs of increasing radii. The DOS 
resonances are seen to be much sharper for metallic CNs 
(m — 3q, q=l,2, ...) than for semiconducting ones (m^ 
3q) in agreement with the fact that this frequency range 
is dominated by the classical Drude-type conductivity 
which is larger in metallic CNs [45| . Figure [4] (b) shows 
the normalized two-particle local photonic DOS functions 
^ taken at the resonance frequencies of the respective 
^'s [Fig. [31(a)], as functions of the distance between the 
two atoms on the axes of the (10,0), (11,0) and (12,0) 
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34 



Interatomic distance (R cn ) 



FIG. 4: (Color online) (a) Transverse local photonic DOS 
functions ^ x (k) for the two- level atom in the centers of the 
four 'zigzag' nanotubes. (b) Normalized two-particle local 
photonic DOS functions £ + (solid lines) and £ _ (dashed lines) 
taken at the peak frequencies of ^(x) [see (a)], as functions 
of the distances between the two atoms on the axes of the 
(10,0) (lines 1; a; = 0.29), (11,0) (lines 2; a; = 0.25) and (12,0) 
(lines 3; a; = 0.24) CNs. (c) Same as in (b) for the two atoms 
on the axis of the (9,0) CN (a; = 0.32). (d) Two-particle DOS 
functions £ + (x) (solid lines) and £~(x) (dotted lines) for the 
two atoms located in the center of the (9,0) CN and separated 
from each other by the distances of 2.1R cn ~7.4A (lines 1), 
i.2R cn » 14.8 A (lines 2) and Q.3R cn «s 22.2 A (lines 3) [shown 
by the vertical lines in (c)]; the dashed line shows £ x (s;) for 
the atom in the center of the (9,0) CN. 



CNs. The values of £ + and £~ are seen to be substantially 
different from each other before they reach their limit 
values of £ + = £~ = £ _L + £" « 2£ x at sufficiently large 
interatomic separations. For metallic CNs the functions 
£ + and £~ exhibit resonator-like behavior, i.e. they vary 
periodically in antiphase almost without damping with 
increasing interatomic separation. In Fig. 3] (c) is shown 
the dependence of ^(x = 0.32) [peak position of £ x in 
Fig. 2] (a)] on the interatomic separation for the atoms 
in the center of the metallic (9,0) CN. As the separation 
increases, £ + and £ _ may differ greatly one from another 
in a periodic way. The maximal difference is |£ + — £~| « 
2(£ ± + f H) « 4£ x - I0 8 , making the mixing coefficients 
C± (t) different and thus resulting in a substantial degree 
of the entanglement of the two spatially separated atomic 
qubits. 

The ensuing spontaneous decay dynamics and atomic 
entanglement are presented in Fig. [5] (a),(b),(c) for the 
atoms separated by the three different distances [shown 
by the vertical dashed lines in Fig. [4] (c)] in the center 
of the (9,0) CN. Atom A is assumed to be initially ex- 
cited while atom B is in its ground state. The entangle- 
ment is seen to reach the amount of 0.5 and to vary with 




FIG. 5: (Color online) (a),(b),(c) Upper decay of the ini- 
tially excited atom A (lines 1) and initially unexcited atom B 
(lines 2), and the entanglement (lines 3), as functions of di- 
mensionless time for the two atoms located in the center of the 
(9,0) CN and separated from each other by the distances of 
2.lR cn « 7.4 A, 4.2_R cn « 14.8 A and 6.3i? c „ » 22.2 A, respec- 
tively [the situation shown in terms of local photonic DOS 
in Figs. |4Tc).(d)]. (d) Short-time and long-time entangle- 
ment evolution of the initially fully entangled atoms sepa- 
rated by 7.4 A; lines 1 and 2 are for C+(0) = 1, C_ (0) =0 and 
C+(0) = 0, C- (0) = 1, respectively. 

time periodically without damping at least for the (rea- 
sonably long) times we restricted ourselves in our com- 
putations. As the interatomic separation increases, so 
the period of the entanglement oscillations does while no 
change occurs in the maximal entanglement. Figure[5](d) 
shows the short-time and long-time entanglement evolu- 
tion when both of the atoms are initially maximally en- 
tangled [C±(0) = 1 while C T (0) = 0] and separated by 
the distance of 7.4 A. The entanglement is slightly larger 
in the case where C+ (0) = 1, C_ (0) = and no damping 
occurs as before. Note that the atoms can be separated 
by longer distances with roughly the same entanglement 
due to the £^ periodicity with interatomic distance. 

4. Exciton-Surface-Plasmon Coupling in Pristine 
Semiconducting CNs 

The total Hamiltonian of an exciton interacting with 
vacuum-type electromagnetic fluctuations (no external 
electromagnetic field is assumed to be applied) on the 
surface of a pristine single- walled semiconducting carbon 
nanotube is of the form (34| 

H = Hp + H ex + H in t , (17) 

where the three terms represent the free field, the free 
exciton, and their interaction, respectively. Here, the 
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FIG. 6: (Color online) Schematic of the transversely quan- 
tized azimuthal electron- hole subbands (left), and the first- 
interband ground- internal-state exciton energy (right) in a 
small-diameter semiconducting carbon nanotube. See text 
for notations. 



second quantized electromagnetic field Hamiltonian is 

poo 

H F = J2 duhujp(n,(j)f(n,u]) (18) 

n J ° 

with the scalar bosonic field operators /'(n, a;) and 
/(n, (J) creating and annihilating, respectively, the 
(vacuum- type) medium- assisted surface electromagnetic 
excitation of frequency u> at an arbitrary point n = R„ = 
{Rcn, <Pn, z n } associated with a carbon atom (represent- 
ing a lattice site) on the surface of the CN of radius Rcn- 
The summation is made over all the carbon atoms, and 
may be replaced by the integration over the entire nan- 
otube surface according to the rule 



i 

So 



1 



2tt 



dip n R C N I dz n 



where S = (i^/^)b 2 is the area of an elementary equi- 
lateral triangle selected around each carbon atom in a 
way to cover the entire surface of the nanotube, with 
6=1.42 A being the carbon-carbon interatomic distance. 

The second quantized Hamiltonian of the free exciton 
(see, e.g., Ref. [77j) i s °f the following form 

H ex = J2 Ef(n)Bl +mJ B mJ =Y,E f (k)Bl f B kJ , 

k,/ 

(19) 

where the operators B a ^ and B n j create and annihilate, 
respectively, an exciton with energy Ef(n) at the lattice 
site n of the CN surface. The index / 0) refers to 
the internal degrees of freedom of the exciton. Alterna- 
tively, Blj = EnBlje 4 '"/^ creates (N is the num- 
ber of the lattice sites) and -Bk./ = ^ annihilates 
the /-internal-state exciton with the quasi-momentum 



k = {A; ¥ ,, k z }, where k v is quantized due to the transverse 
confinement effect and k z is continuous. The exciton to- 
tal energy is given by 



E f (k) = Egl(k v ) 



h 2 k 2 z 

2M PX 



(20) 



with the first term representing the excitation energy, 



EJK 



( f) 

E^ ' , of the /-internal-state ex- 

(/) 



citon with the (negative) binding energy E b , created 
via the interband transition with the band gap E g (k ip ) = 
£e(k v ) + £h(k v ), where e e .h are transversely quantized 
azimuthal electron-hole subbands (see the schematic in 
Fig. [SI). The second term represents the kinetic energy 
of the translational longitudinal movement of the exci- 
ton with the effective mass M ex = m e + nth , where m e 
and rrih are the electron and hole effective masses, re- 
spectively. The two equivalent free-exciton Hamiltonian 
representations are related to one another via the obvi- 
ous orthogonality relationships ^ n e - l ( k ~ k )' n /jV = 
and e ~ l ( n ~ m ) ' k /iV = (5 nm with the k-summation run- 
ning over the first Brillouin zone of the nanotube. The 
bosonic field operators in Eq. (fT5)) are transformed to the 
k-representation in the same way. 

The most general (non-relativistic, electric dipole) 
exciton-electromagnetic-field interaction on the nanotube 
surface can be written in the form 34 , 48 , 5(| 



H;, 



H 



(i) 



H 



(2) 



(21) 



m e c 



E A ( n )' 



Pn - ^A(n 



^d n - V n 0(n) . 



where p n = 5Z/(0|pn|/)5 n ,/ + h.c. is the total electron 
momentum operator at the lattice site n under the optical 
dipole transition resulting in the exciton formation at the 
same site, d n = J^(0|d n |/)i? n ./ + h.c. is the correspond- 
ing transition dipole moment operator [related to p n via 
the equation (0|p n |/) = im e Ef(n)(0\d n \f) /he], c and e 
are the speed of light and the electron charge, respec- 
tively. The vector potential operator A(n) (the Coulomb 
gauge is assumed) and the scalar potential operator ip(n) 
represent, respectively, the nanotube's transversely po- 
larized surface electromagnetic modes and longitudinally 
polarized surface electromagnetic modes which the exci- 
ton interacts with, 



c * -L 

A(n) = / du — E (n, u>) + h.c. 

o * w 



■^TrhioRea zz (RcN, w) ± G zz {n, m, w)/(m, u) + h.c, 



V n ip(n) = du>~E (n, ui) + h.c. 



|>OC 



4iuj 



y / irhujRc<j zz (R C n,u) n G zz (n,m,uj)f(m,uj) + h.c., 
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with the total electric field operator given by 



E(n): 



dw 



E (n,w) 



•E"(n,w) 



h.c. 



and ± ^G zz (n, m, oj) representing the zz-component of 
the transverse (longitudinal) Green tensor of the electro- 
magnetic subsystem (the nanotube). This Green tensor 
is different from that used for atomically doped CNs in 
Sections 2 and 3 above since now there is no additional 
atomic source in the problem. The Green tensor is de- 
rived by expanding the solution of the Green equation in 
cylindrical coordinates and determining the Wronskian 
normalization constant from the appropriately chosen 
boundary conditions on the CN surface 0, [Hol.]7§| . 

Equations (fT7|) - (f2"Tj) form the complete set of equations 
describing the exciton-photon coupled system on the CN 
surface in terms of the electromagnetic field Green tensor 
and the CN surface axial conductivity. 

It is important to realize that the transversely polar- 
ized surface electromagnetic mode contribution to the in- 
teraction Hamiltonian from Eq. pTj) (first term) is neg- 
ligible compared to the longitudinally polarized surface 
electromagnetic mode contribution (second term). The 
point is that, because of the nanotube quasi-one- 
dimensionality, the exciton quasi-momentum vector and 
all the relevant vectorial matrix elements of the momen- 
tum and dipole moment operators are directed predom- 
inantly along the CN axis (the longitudinal exciton). 
This prevents the exciton from the electric dipole cou- 
pling to transversely polarized surface electromagnetic 
modes as they propagate predominantly along the CN 
axis with their electric vectors orthogonal to the prop- 
agation direction. The longitudinally polarized surface 
electromagnetic modes are generated by the electronic 
Coulomb potential (see, e.g., Ref. [ZSj]), and therefore 
represent the CN surface plasmon excitations. These have 
their electric vectors directed along the propagation di- 
rection. They do couple to the longitudinal excitons on 
the CN surface. Such modes were observed in Ref. [56| . 
They occur in CNs both at high energies (well-known n- 
plasmon at ~ 6 eV) and at comparatively low energies of 
~ 0.5 — 2 eV. The latter ones are related to transversely 
quantized interband (inter-van Hove) electronic transi- 
tions. These weakly-dispersive modes [5^, H(| are simi- 
lar to the intersubband plasmons in quantum wells [8lj |. 
They occur in the same energy range of ~ 1 eV where the 
exciton excitation energies are located in small-diameter 
(< 1 nm) semiconducting CNs [5?], HU . 

The Hamiltonian (fT7| - ((2T|) can be diagonalized exactly 
using Bogoliubov's canonical transformation technique 
(see, e.g., Ref. [52]), to yield 34 1 



H= J2 ^(k)H(k)Uk)+E 0l (22) 

k, fj,=l,2 



4(k) = J2 [u*(k, Uf )B Kf -iv(k, W/ )B! 



iM 1 



-k,/ 



t (k,o;)/(k,a;)-<(k,a;)/t(-k,a;) 



annihilates and £t(k) = [£ M (k)]t creates the coupled 
exciton-plasmon excitation of branch /i (= 1,2), and 
are the (appropriately chosen) canonical transforma- 
tion coefficients, ujf — Ef/Ti. The "vacuum" energy Eq 
represents the state with no exciton-plasmons excited in 
the system, and ?kJ M (k) is the exciton-plasmon energy 
given by the solution of the following (dimensionless) dis- 
persion relation 34, 3f| 



EL 
2tt 



dx x f q (x) ■ 



p{x) 



= 0. 



(23) 



Here, 



the function 



fc^(k) 
270 



270 



4|4| 2 x 3 
3?ic 3 



2to 

h 



represents the (dimensionless) exciton spontaneous decay 
rate, where d{ =Y^ a {Q\{dn)z\f) is the longitudinal exci- 
ton transition dipole moment matrix element, and the 
function 



p(x) 



3S 



Re- 



lQiraR 2 CN cr zz (x) 



(24) 



stands for the surface plasmon DOS responsible for the 
exciton decay rate variation due to its coupling to the sur- 
face plasmon modes. In Eq. (f24l) . a — e 2 /he— 1/137 is the 
fine-structure constant and a zz = 2 r r:fia zz /e 2 is the dimen- 
sionless CN surface axial conductivity per unit length. 
Note that the conductivity factor in Eq. (|24[) equals 



Re- 



1 



4ac / h 



Im- 



1 



where the new operator 



cr zz (x) R C n\2iqxJ e zz (x)-l 

in view of the equation 

_ iu)(e zz -l) 
AttSpt 

representing the Drude relation for CNs, where e zz is the 
longitudinal dielectric function, S and pr are the surface 
area of the tubule and the number of tubules per unit vol- 
ume, respectively 0, Bill Is3- This relates very closely 
the surface plasmon DOS function (|24|) to the loss func- 
tion — Im(l/e) measured in Electron Energy Loss Spec- 
troscopy (EELS) experiments to determine the properties 
of collective electronic excitations in solids [56fl . 
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FIG. 7: (Color online) Fragments of the real/imaginary con- 
ductivity and the surface plasmon DOS (a), and the lowest 
bright exciton dispersion (b) when coupled to plasmons in the 
(11,0) CN. Dimensionless energy is defined as [Energy] /2 r yo- 
See text for dimensionless quasi-momentum. 



tion]. The peak of the DOS function is seen to exactly 
coincide in energy with the zero of Im[<7 2Z (a;)] {or the zero 
of Re[e zz (a;)]}, clearly indicating the plasmonic nature of 



the CN surface excitations under consideration 
Also seen in Fig. that the interband plasmon excita- 
tions occur in CNs slightly above the first bright exciton 
excitation energy [70j. This is a unique feature of the 
complex dielectric response function — the consequence 
of the general Kramers-Kronig relation [62| . 

The sharp peak structure of the surface plasmon DOS 
p(x) may be used to solve the dispersion equation (j23|) 
analytically. Using the Lorentzian approximation 



p(x) 



p{x p )Ax 2 v 



x p ) 2 + Axj 



with x p and Ax p being, respectively, the position and the 
half-width-at-half-maximum of the plasmon resonance 
closest to the lowest bright exciton excitation energy, the 
integral in Eq. (f23|) simplifies to the form 



dx 



xf f (x)p(x) F(x p )Ax 2 p 



dx 



F(x p ) 



arctan 



V Ax , 



(x — x p ) 2 + Ax 2 



where F(x p ) — 2x p T( ) {x p )p{x p )/-K. This is valid for all 
x^ in Eq. (|23p apart from those located in the narrow in- 
terval {x p — Ax p , x p + Axp) in the vicinity of the plasmon 
resonance, provided that the resonance is sharp enough. 
Then, the dispersion equation becomes the biquadratic 
equation for x^ with the following two positive solutions 
of interest (the dispersion curves) 



Figure [2] above shows the low-energy behavior of the 
functions a zz (x) and Re[l/er zz (x)] for the (11,0) nan- 
otube taken as an example (Rcn — 0.43 nm). The func- 
tion Re[l/(7 zz (:r)] is only non-zero when the two condi- 
tions, lm[a zz (x)]_= and Re[a zz (a;)] — > 0, are fulfilled 
simultaneously [73|, 0, 81 1. These result in the peak 
structure of Re(l/tr zz ) and of the surface plasmon DOS 
function (|2"4")l . respectively. The plasmon DOS peaks rep- 



1 



resent the surface plasmon modes associated with the 
transversely quantized interband electronic transitions in 
CNs [80|| (transitions in between the azimuthal electron- 
hole subbands in the left panel of Fig. [5]). 

Figure (a) shows the lowest-energy plasmon DOS 
resonance of the (11,0) CN [H as given by Eq. ([23]). 
The corresponding fragments of the nanotube's real and 
imaginary conductivities are also shown there for com- 
parison. The lower dimensionless energy limit is set up 
to be equal to the lowest bright exciton excitation energy 
[E exc = 1.21 eV (x = 0.224) for the (11,0) CN as reported 
in Ref. 57] by directly solving the Bethe-Salpeter equa- 



±- J{e)-xlf + Fpe s 



(25) 



Here, F p 



4F(x p )Axp(ir — Axp/xp) (the arctan-function 
is expanded into series to linear terms in Ax p /x p < 1). 

The two dispersion curves of Eq. (|25l) are shown in 
Fig. [7] (b) for the (11,0) nanotube [35| as functions 
of the dimensionless longitudinal quasi-momentum. In 
these calculations, the dipole interband transition ma- 
trix element in Tq(x p ) was estimated from the equation 
\d{\ 2 = 3HX 3 1 4:T™ d (according to Hanamura's general 
theory of the exciton radiative decay in spatially con- 
fined systems [83]), where T™ d is the exciton intrinsic 
radiative lifetime, and A = lixchjE with E being the ex- 
citon total energy given by Eq. (f2"0"l) . For zigzag-type nan- 
otubes, the first Brillouin zone of the longitudinal quasi- 
momentum is given by —2nh/3b < Hk z < 2Trh/3b [3,0]. 
The total energy of the ground-internal-state exciton can 
then be written as E = E exc + (27r?i/36) 2 i 2 /2M ex with 
— 1 < t < 1 representing the dimensionless longitudinal 
quasi-momentum. The lowest bright exciton parameters 
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r rad 

ex 



14.3 ps, and M ex = 0.44m (mo 
is the free-electron mass) were used for the (11,0) CN as 
reported in Ref.[13] from the numerical solution of the 
Bethe-Salpeter equation. 

The dispersion curves in Fig. [7] (b) demonstrate a 
clear anticrossing behavior with the (Rabi) energy split- 
ting ~ 0.1 cV. This indicates the formation of the 
strongly coupled surface plasmon-exciton excitations. 
Note that here the strong exciton-plasmon interaction 
is supported by an individual quasi-lD nanostructure, 
a single-walled (small-diameter) semiconducting carbon 
nanotube, as opposed to the artificially fabricated metal- 
semiconductor nanostructures studied previously 
l6ll | where the metallic component normally carries the 
plasmon and the semiconducting one carries the exciton. 
It is also important that the effect comes not only from 
the height but also from the width of the plasmon res- 
onance as is seen from the definition of the F p factor in 



Eq. (|25p . In other words, as long as the plasmon reso- 
nance is sharp enough (which is always the case for in- 
terband plasmons), the effect is determined by the area 
under the plasmon peak in the DOS function (|24p rather 
then by the peak height as one would expect. 

The theory of Section 2 can be applied to obtain the 
optical absorption lineshape (optical response) for the 
surface exciton-plasmons in CNs. A slight difference now 
is that the exciton-phonon scattering should also be taken 
into account as a concurrent process to scatter the exci- 
ton in addition to the exciton-plasmon scattering. In the 
relaxation time approximation for the exciton-phonon 
scattering, the (dimensionless) exciton absorption line- 
shape function I(x) in the vicinity of the plasmon reso- 
nance is then of the form [35j 



J (e/) [(x - e f f + Ax 2 p 



[(x-Sff-Xy^ + (x-e f r(Axl+As}) 



(26) 



I(x) = 



with 



Io(e f ) = ^n(e f )p(Ef), X f = ^4irAx p Io(s f ), 

and Ae f — hjl^T^ being the exciton energy broadening 
due to the exciton-phonon scattering with the relaxation 
time T ph - 

Figure [8] shows the calculated exciton absorption line- 
shapes as the exciton energy is tuned to the nearest plas- 
mon resonance of the (11,0) CN. In these calculations, the 
phonon relaxation time r p h = 30 fs was used as reported 
in Ref. [3l|. The line (Rabi) splitting effect is seen to be 
~0.1 eV, indicating the strong exciton-plasmon coupling 
with the formation of the mixed surface plasmon-exciton 
excitations. The splitting is not masked by the exciton- 
phonon scattering, and is expected to be larger in smaller 
diameter nanotubes. 

Obviously, the formation of the strongly coupled mixed 
exciton-plasmon states is only possible if the exciton to- 
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FIG. 8: (Color online) Exciton absorption lineshape as the 
exciton energy is tuned to the nearest plasmon resonance 
[vertical dashed line here; see Fig. [7] (a)] in the (11,0) CN. 
Dimensionless energy is defined as [Energy] /2"/o- 



tal energy is in resonance with the energy of an inter- 
band surface plasmon mode. The exciton energy might be 
tuned to the nearest plasmon resonance in ways used for 
the excitons in semiconductor quantum microcavities — 
thermally (by elevating sample temperature) [H^, 54 , 5f| , 
and/or electrostatically [H,l83,[8l,|8§| (via the quantum 



confined Stark effect with an external electrostatic field 
applied perpendicular to the CN axis). The two possi- 
bilities influence the different degrees of freedom of the 
quasi-lD exciton — the (longitudinal) kinetic energy and 
the excitation energy, respectively [see Eq. ([20])]. In the 
latter case, in spite of the fact that the cylindrical sur- 
face symmetry of the excitonic states brings new pecu- 
liarities, the general qualitative behavior of the quantum 
confined Stark effect in CNs should be similar to what 
was previously observed and theoretically analyzed for 
semiconductor quantum wells [86* 



87] 



One should ex- 
pect that the exciton excitation energy E exc and the in- 
terband plasmon energy E p both shift to the red due to 
the decrease in the CN band gap E g as the perpendicu- 
lar electrostatic field increases. However, the exciton red 
shift is expected to be much less due to the decrease in 
the absolute value of the exciton (negative) bind- 
ing energy (9p| . which is estimated to be ~ 0.5E^ XC in 
small-diameter carbon nanotubes [2^, [2l|, [22|, [25| and 
thus contributes largely to the exciton excitation energy. 
So, the exciton excitation energy E exc and the interband 
plasmon energy E p will be approaching as the field in- 
creases, bringing the total exciton energy [see Eq. (|20|1 ] 
in resonance with the plasmon mode due to the non-zero 
longitudinal kinetic energy term at finite temperature. 
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FIG. 9: (Color online) Calculated dependences of the first 
bright exciton parameters in the (11,0) CN on the electro- 
static field applied perpendicular to the CN axis. Dimension- 
less energy is defined as [Energy]/ '270. See text for notations. 



Figure [9] shows the results of the recent calculations 
of the quantum confined Stark effect for the (11,0) nan- 
otube [90J, which confirm the qualitative expectations 
just discussed. More details on these calculations and the 
complete theory of the quantum confined Stark effect in 
CNs will be published elsewhere. 

5. Conclusion 

Similar to semiconductor microcavities and photonic 
band-gap materials, carbon nanotubes may qualitatively 
change the character of the atom-electromagnetic-field 
interaction, yielding the strong atom-field coupling. For 
atoms (ions) doped into small-diameter (~ 1 nm) nan- 
otubes, the effect of the optical absorption line splitting 
(Rabi splitting) of ~ 1 eV is predicted in the vicinity 
of the energy of the atomic transition. This is at least 
three orders of magnitude larger than the Rabi splitting 
of ~ 140—400 /xeV observed recently in the photolumines- 
cence experiments for exciton-polaritons in quantum dots 
in semiconductor microcavities 53|, 54 , [Hi] . The larger 
effect in atomically doped CNs comes from the facts that 
(i) typical atomic binding energies are at least three or- 
ders of magnitude larger than typical excitonic binding 
energies in solids, and (ii) the effective mode volume of 
the medium-assisted vacuum-type electromagnetic field 
in the CNs is small enough to achieve Rabi splitting of 
~ 1 eV in 1 nm-diameter nanotubes. 

Rear-earth Eu 3+ complexes in hydrocarbon solvents 
might be good candidates to test this prediction. They 
are known to be excellent probes to study quantum optics 
phenomena in spatially confined systems such as dielec- 



tric microspheres 84| , or photonic crystals [85| , owing to 
the dominant narrow, easily detectible Dq—* 7 F<2 dipole 
transition between two deep-lying electronic levels of eu- 
ropium that essentially create a two-level system. 

Metallic carbon nanotubes of ~ 1 nm diameter are 
demonstrated theoretically to be very efficient, even at 
room temperatures, to entangle a pair of the spatially 
separated atomic qubits located in the middle of the nan- 
otube, provided they interact with the same resonance 
of nanotube's local photonic DOS. This latter condition 
is automatically satisfied for identical atoms doped into 
the nanotube. The entanglement greatly exceeds 0.35 
(known to be the maximal value one can achieve in the 
weak atom- field coupling regime 65]), and persists with 
no damping for very long times. Such a behavior is 
attributed to the strong atom-field coupling and elec- 
tronic structure resulting in the resonator-like distance 
dependence of the two-particle local photonic DOS in 
metallic CNs. Thus, similar to quasi-OD excitonic po- 
laritons in quantum dots in semiconductor microcavities 
suggested recently to be a possible way to produce the 
excitonic qubit entanglement [43j |. small-diameter metal- 
lic atomically doped carbon nanotubes might represent 
a novel approach towards solid-state quantum informa- 
tion transfer over long distances — centimeter-long dis- 
tances, as a matter of fact, since centimeter-long small- 
diameter single- walled CNs are currently technologically 
available 0, 0]. 

In the current version of the theory outlined here, all 
the possible decoherence mechanisms in the two-qubit 
system coupled to nanotube's surface electromagnetic 
mode are taken into account in terms of the relaxation 
time approximation in the CN surface conductivity. The 
conductivity appears in the photonic DOS and plays a 
crucial role in the entanglement of formation. It is also 
important that the method of entanglement preparation 
described here does not involve the spin degrees of free- 
dom — just a normal, but strong, electric dipole coupling 
to the same local (single-particle) photonic DOS reso- 
nance of the nanotube is required to entangle the pair of 
spatially separated atomic qubits, as opposed to the spin 
entanglement studied recently in other fullerene-based 
models 42j. The outlined entanglement scheme can be 
further generalized to the multi-atom entanglement via 
the nearest neighbor pairwise quantum correlations, thus 
challenging novel applications of atomically doped CNs 
in quantum information science. 

In pristine small-diameter (< 1 nm) semiconducting 
CNs, the strong exciton-surface-plasmon coupling effect 
is demonstrated with the characteristic exciton absorp- 
tion line (Rabi) splitting ~ 0.1 eV. This is almost as 
large as the typical exciton binding energies in such CNs 
(~ 0.3-0.8 eV [H El HI Hi|), and of the same or- 



der of magnitude as the exciton-plasmon Rabi splitting 
in organic semiconductors (~ 180 meV [59(). Also, this 
is much larger than the exciton-polariton Rabi splitting 
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in semiconductor microcavities (~ 140 — 400 fj,eV 5^, 54, 
55J), or the exciton-plasmon Rabi splitting in hybrid 
semiconductor-metal nanoparticle molecules [6(|. Also 
important is that the strong exciton-plasmon coupling 
effect takes place in an individual carbon nanotube as 
opposed to artificially fabricated hybrid plasmonic nanos- 
tructures just mentioned. To bring the exciton in reso- 
nance with nanotube's nearest plasmon mode, the exci- 
ton energy can be tuned in ways used for the excitons in 
semiconductor quantum microcavities thermally (by 



elevating sample temperature) [53l . |54| . |55| . and/or elec- 
trostatically 0, [St], lH, [89| (via the quantum confined 
Stark effect with an external electrostatic field applied 
perpendicular to the CN axis). The two possibilities af- 
fect the different degrees of freedom of the quasi-lD exci- 
tons in CNs — the (longitudinal) kinetic energy and the 
excitation energy, respectively. 

To conclude, the surface electromagnetic field phenom- 
ena reviewed here have a great potential to open up new 
paths for the development of the CN based tunable opto- 
electronic device applications in areas such as nanopho- 
tonics, nanoplasmonics, cavity quantum electrodynam- 
ics, and quantum information science. 
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